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Abstract 


The problem of heat conduction on networks of multiply connected rods is solved 
by providing an explicit solution of the one-dimensional heat equation in each domain. 
The size and connectivity of the rods is known, but neither temperature nor heat flux 
are prescribed at the interface. Instead, the physical assumptions of continuity at the 
interfaces are the only conditions imposed. This work generalizes that of Deconinck, 
Pelloni, and Sheils [S] for heat conduction on a series of one-dimensional rods connected 
end-to-end to the case of general configurations. 


1 Introduction 


Suppose n is a set of vertices, and A is a set of directed edges between those vertices so 


that {V, E) is a finite connected directed graph. Associated with each edge r G A is the 
length Lr G (0, oo]. Suppose additionally if L,, = oo, then the vertex at which r terminates 
has no other edges. For each edge r £ E, we dehne the open interval Clr = (0,^^). We 
consider the graph to represent a configuration of narrow metal rods joined at the vertices 
in perfect thermal contact. In this paper we study the distribution of heat on such a 
configuration. 

In order to derive the proper equations and interface conditions, we generalize what is 
outlined in [T2| for problems with multiple rods. Assume each edge r G A is a rod made 
of some heat-conducting material with density pr{x) and unit cross-sectional area. We 
assume the surface of the rod is perfectly insulated so no heat is lost or gained through 
this surface. If we consider an inhnitesimal section of length dx for x £ then dC"", 
the heat content in the section, is proportional to the mass and the temperature, q^{x,t). 
That is 


dC"’(t) = Cr{x)pr{x)q^{x, t) dx. 


where Cr-(x) is the specific heat on r and x £ Q.r- Thus, the total heat content in the 
interval xi ^ x ^ X2 where xi,X2 G Qr is 
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Fourier’s Law for heat conduction im states that the rate of heat flowing into a body 
is proportional to the area of that element and to the outward normal derivative of the 
temperature at that location. The constant of proportionality is the thermal conductivity, 
kr{x). In our example, the net inflow of heat through the boundaries xi and X2 is 


R{t) = kr{x2)ql{x2,t) - kr{xi)ql{xi,t), (1.1) 

where xi,X 2 G Hr- Conservation of heat implies = R{t) along each rod. That is, 

( 1 . 2 ) 


— J Cr{x)pr{x)q'~{x,t)dx = kr{x 2 )qlix 2 ,t) - kr{xi)ql{xi,t), 


for every r. This is a typical conservation law. On each rod we assume constant material 
properties. That is Cr{x) = Cr,Prix) = pr, and kr{x) = kr for x G 11^ and 0 < t < T. 


Expressing R{t) as the integral of a derivative and rewriting (1.2), we have 
j (^CrPrqlix,t)dx --^{krqlixR))^ dx = 0, 


(1.3) 


for xi,X 2 G Hr- Since this is true for any xi and X 2 in this range, it follows that the 
integrand must vanish. That is, 

- —ql^{xR) = 0 . 

Cj' P'f' 

Next, we scale x such that x = for x G Hr- Note that such a scaling also affects 

Crpr 

Lr and and in what follows we assume all quantities are properly scaled. Dropping 
the ", and defining cr^ = krCrPr as the (scaled) thermal diffusivity in rod r we have 


r _ ^2 r 

^rHxx') 


g^(x,0) = gS(x), 


X G Dr-) t > 0 , 
X G Dr-) 


r £ E, (1.4a) 

r £ E, (l-4b) 


where q^ £ 5 (D,.) is given initial data for each r £ E with S (D) a Schwartz space of 
smooth, rapidly decaying functions restricted to the closure of the interval D. Wherever 
D is finite, S (D) = C°° (D). 

We also assume that the rods at each vertex v are in perfect thermal contact [3]. The 
temperature at x = 0 of the rods emanating from v, q^{0, t), is the same as the temperature 
at X = Lr of the rods terminating at v, q'^{Lr,t). That is. 


3 0^ : [0, T] —)> M : V r G E, if r emanates from v then q'^{0,t) = 6r{t), t £ [0,r], 


(1.5) 


y r £ E, if r terminates at v then q^{Lr,t) = 9v{t), t £ [0,T]. 

Note that the functions 9^ are not data of the problem; they appear in the interface 
condition only for notational convenience. If vertex v has p edges then (1.5) gives p — 1 
conditions. 

In order to hnd a condition on the set of spatial derivatives of q^ we require an appro¬ 
priately scaled (1.2) valid on a region centered at the vertex v. Consider a ball centered 
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at a vertex v, with radius e > 0 sufficiently small so that all rods are straight lines within 
the ball and no other rods intersect the ball. Then, similar to Equation (1.1) we have 

R{t)= ^ULr-e)qliLr-e,t). ( 1 . 6 ) 


sGE: 

s emanates 
from V 


r£E\ 

r terminates 
at V 


By conservation of energy, 

/ 


d 

dt 


\ 


c^Pslim I g®(x,t)dx+ lim I^ q^{x,t)dx 

sGE: 
emanates 


0 


from V 


r£E: 

r terminates 
at V 


Lr — £ 


^s9x(e,t)-hm V ajqliLr - e,t). 


seE: 
s emanates 
from V 


e^O 


reE: 
r terminates 
at V 


The left-hand-side of ( |1.7[ ) is zero by ( |1.5[ ). This implies 

Y Y ^hl{Lr,t) = 0. 


sGE: 

s emanates 
from V 


reE: 

r terminates 
at V 


(1.7) 


( 1 . 8 ) 


Thus the appropriate sum of the heat flux is continuous across the interface. 

For each finite endpoint vertex v £ V which has only one finite-length edge £ E 
(and no infinite-length edges) we prescribe a general Robin boundary condition 


/3-oq^^{X,t)+/3ld^q'^-{X,t) = Mt), t £ [0,r], (1.9) 


where G C'°°[0, T] is known boundary datum and (3^,(3^ £ M where X = 0 if 
emanates from v, and X = if terminates at u. If the graph has rrif finite edges 
and nii infinite edges, then (1.5), (1.8), and (1.9) together prescribe a total of 2m/ -|-m* 
boundary and interface conditions. 

We insist that the initial data are mutually compatible according to the interface 
conditions. That is, we require for each interface vertex v, 


Y ^sdxqoi^)- Y cr'^d^qoiLr) = 0, 

sGE: vGE: 

s emanates r terminates 

from V at i; 

V s G E, if s emanates from v then (?o(0) = ^i;(0), 
y r £ E, if r terminates at v then q^iLr) = 0^(0). 


We also require the initial and boundary conditions are compatible in the sense that for 
each finite endpoint vertex v with edge r^. 


/3oX''W+/3r5x<?S''W = ^(o)- 


Deconinck, Pelloni, and Sheils [5] apply the Fokas method, alternatively called the 
unified transform method, to a subclass of the problems described above in which each 
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vertex has at most two edges: a configuration of metal rods with different diffusivities 
joined end-to-end along a line. They derive the so-called global relation for each rod using 
Green’s theorem applied to the space-time domain, and follow the approach of [7j to obtain 
an implicit integral representation for the solution, known as the Ehrenpreis form [H]. 

This representation depends upon (time-transforms of) all boundary and interface 
values diq'^{0,t), diq'^{Lr,t) for j G {0,1}, r € E, t ^ [0,T], so cannot be considered an 
effective solution representation. Using the boundary and interface conditions together 
with the global relation, the authors are able to formulate and solve a system of linear 
equations for these transformed boundary values in a number of particular cases: two 
or three semi-infinite, or finite rods. It is clear that the method used would work for a 
configuration with more than three rods, but the linear algebra becomes cumbersome to 
solve by hand. 

The principal contribution of the present work is to show that the Fokas method 
may also be applied to the more general configurations of rods described above. The 
heat equation on graphs is important in engineering applications such as building and 
analyzing heat exchanger networks [ISlEs] and “channelling” heat through materials of 
differing conductivity [3] . 

1.1 Organization of paper 

In Section]^ we derive a global relation on each domain Qr x (0,T) and use it to obtain 
in implicit integral representation for the solution in Ehrenpreis form. 

In the proceeding sections we present explicit solutions for several cases of relatively 
simple configurations of rods. Specifically, in Sectionwe study the problem of arbitrarily 
many semi-infinite rods joined at a single vertex. In Section we consider arbitrarily 
many finite parallel rods between a pair of common vertices, and in Section we analyze 
arbitrarily many finite rods emanating from a single common vertex, and present an 
example of such a network of three rods. 

We conclude with some remarks on the applicability of the method presented herein 
to related problems. 

2 Implicit integral representation of the solution 

2.1 Global relation 

The so-called global relation is an essential tool in the Fokas method. As our problem 
is posed on domains fl,. for each r € E we will derive \E\ global relations. Each global 
relation links the boundary and interface values with the initial value for a particular rod, 
and holds for all values of the complex spectral parameter A within a half or full plane. 
In contrast with [5j, but following the innovation of mm, we scale A in a way that will 
simplify later calculations. 

For {x,t) G Qj. and A G C, we dehne the functions 

Xr{x,t,X) = 

Yr{x,t,X) = ^(Jrix + cr^dx) q''{x,t), 
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which have the property {dtXr — dxYr){x, t,X) = 0 if and only if q'^ satisfies equation 
Suppose that satisfies the partial differential equation and the initial condition 
Then Green’s theorem [T] applied to Qrit) = (0,Lr.) x (0,t) yields 

0 = / (9tX — s, A) dx ds = f (Ydt + Xdx). 

jQrit) J 


Rearranging equation (2.1) and using the notation 


foW = J 

f e ^^^qo{x) dx, 

0 

AeC 

fiX;t) = J 

1 e““ q{x,t)dx, 

0 

ft 

AeC 

55(w,t) = J 

1 e‘^^q^iO,s)ds, 

0 

ft 

w E C 

gl{u,t) = J 

1 e‘^^ql{0,s)ds, 

0 

ft 

w E C 

II 

1 e‘^’^q'^{Lr,s) ds, 

0 

ft 

cu E C 

II 

1 e‘^^ql{Lr,s)ds, 

0 

oj eC 


t E [o,r], 
t E [o,r], 
t E [o,r], 
t E [o,r], 
t E [o,r], 


yields the global relation 




(crriA/io(A^, t) + cr‘^h'[{X^,t)) , 


for A E C and t E [0, T] which holds for each r & E such that Lr < oo. Similarly, 


A 


ari\go{X\t) + a'^gl{X^,t) = , AeC ,tE[0,r], 

for each r € E such that Lr = oo. 

2.2 Implicit integral representation of the solution 


Rearranging the global relation (2.8b), and taking inverse Fourier transforms ir 


and (2.8b) respectively we have 




27r 


/. 


Jx\lar-Yt (.;\„r(\2 


— OO 
OO 


{iXgQ{X^,t)+ argl{X'^,t)) dX 

^iXlar{x-Lr)-XH (iAh[i(A2,t) +U,/l^(A2,t)) dA, 


( Mbt - 

( 2 . 1 ) 

( 2 . 2 ) 

(2.3) 

(2.4) 

(2.5) 

( 2 . 6 ) 
(2.7) 

(2.8a) 

(2.8b) 


(2.8a) 


(2.9a) 
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and 


1 r 




- ^ y + ^rgliX^t)) dA. 


(2.9b) 


We define the sectors 



Re(A) 


Figure 1: The sectors and D , and their oriented boundaries. 


Z)± = {A e C : f < arg(±A) < f} 


( 2 . 10 ) 


as shown in Figure The integrand of the second integral in (2.9b) is entire and decays as 
A —)• oo from within C'^\D~^. Hence, by Jordan’s Lemma [T] we can replace the contour of 
integration of the second integral by fgjj+ ■ Further, by analyticity we make an additional 
contour deformation from dD^ to dD^, where 


= {A E D 


± 


|A| >i?}. 


( 2 . 11 ) 


and i? > 0 is an arbitrary constant. An appropriate (sufficiently large) value of this 
constant may be chosen for any individual problem as in Figure 

We have arrived at the implicit integral representation of the solution known as the 
Ehrenpreis form. For each r € E such that L/j = oo, 

’-^'^"*g5(A) dA 

(2.12a) 

{iXg^,{X^,t) + argl{X^t)) dA. 
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Re(A) 


Figure 2: The sectors and Dj^, and their oriented boundaries. 


Similarly, for each r £ E such that Lr < oo, 
/-oo 






— OO 

_ ^ f ^ixXjar—y^t 

27r JdD+ 

1 r 


e—'"' - “ (^A5^S(A^^) + (Tr5^(A^^)) dA 


(2.12b) 


27r 


'dD- 


^iix-u)xiar-xH (^iXhl^^x\t)+arh[{X^,t)) dA. 


The process of using the initial, boundary, and interface conditions to determine the 
t-transformed boundary values is known as a generalized spectral Dirichlet-to-Neumann 
map; generalized because the map is from the known boundary and initial data to the 
unknown boundary and interface values, rather than from 0**^ order boundary data to 
1®* order boundary values; spectral because the map is formulated in spectral space in¬ 
stead of coordinate space. In each of Sections mu we implement the generalized spectral 
Dirichlet-to-Neumann map for particular classes of configurations of rods. 


3 m semi-infinite rods 

In this section, we consider the case of m G N semi-infinite rods of differing thermal 
diffusivity, joined at a single vertex. The graph is as shown in Figure Of the m -|- 1 
vertices, m are each the terminus of a single edge, and all m edges emanate from the other 
vertex. As \E\ = m, we find it notationally convenient to identify the set of edges E with 
the set {1, 2,..., m}. The single interface vertex is denoted v. 
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Figure 3: m semi-infinite rods. 


In what follows we find explicit representations of the t-transformed interface values 


appearing in Equation (2.12a) in terms of only the initial data and certain remainder 


terms. We then show the remainder terms do not to contribute to 


For each 1 ^ r ^ m, we evaluate the global relation (2.8b) at A i—>■ —A to obtain the 
relations 


- ari\gl{X^,t) + olgl{\^,t) = 


-A 


(7 T 


- 


-A 


(Tip 




(3.1) 


for all A G C^, where we have used that gj{X'^,t) is invariant under A —)• —A. This 
is essential because we require equations for g‘^AX^,t) valid for A G C C^, whereas 


the original global relation equation (2.8b) was valid only for A G C . Equation (3.1) 


represents a system of m linear equations in the 2m unknown, t-transformed interface 
values g'j{X^,t): for 1 ^ r ^ m, j G {0,1}. 


Appropriately transforming the interface conditions (1.5) and (1.8) for this problem 
we find 

go{X^,t) ■= gl{X^A) = f e^^^e^{s)ds, 1 ^ r ^ m, (3.2) 

JQ 

m 

Y,algl{X\t) = Q, (3.3) 

r=l 


for all A G C. Equation (3.2) reduces the number of unknown functions from 2m to m -|-1. 
Equations (3.1), (3.3) provide m -|- 1 linear equations in these unknown functions. We 

AX = Y, (3.4) 


represent this system as a matrix: 
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where 


(3.5) 

(3.6) 


^ = (5o(A^^),o■?5^(A^,^),...,^7^fff(A^,^))"^, 

/ 0 1 1 ••• 1 \ 

—i\ai 1 0 • • • 0 

^ — I —i\a2 0 1 • • • 0 


\-i\am 0 0 ••• 1 / 


_gA2i.m 


T 


(3.7) 


For notational convenience, the entries of Y are denoted 17, where r is enumerated 
0,..., m. Similarly, the rows of the other matrices are counted from 0 to m. 

Clearly, the “data” vector Y depends not only upon the initial data, but also upon 
the Fourier transforms of the solutions at time t which are not given data of the problem. 
Nevertheless, we proceed to solve the system as if Y is composed purely of known data 
and show afterwards that terms involving make no contribution to the solution 

representation. 

Subtracting the sum of rows 1,..., m from row 0 in matrix A it is immediate that 


m 

det^ = i\ (Jr- 

r=\ 


Let Aj be the matrix formed by replacing column j of A with the column vector Y . By 
Cramer’s rule. 


ffo(A^^) 


det ^0 
det^ ’ 


^/^(A^^) 


det Ar 
det A ’ 


1 ^ r ^ m. 


Trivially, 

m 

det^o = — Yp. 

p=i 

To find det Ar we subtract the row for j = 1, 2,..., r — 1, r + 1,..., m from the first 
row. Then, for each j = r, j = r — 1, ... , j = 2, we switch the and {j — 1)®* rows and 
switch the and {j — 1)®* columns. The resulting matrix is 


m 

iX cjp 

m 

0 •• 

■ "1 

P=1 

p=l 



p¥=r 

P¥=r 



'l\CT 

Yr 

0 •• 

• 0 

%y^cT 

Yr> 

1 •• 

• 0 

rpf f 

Yrn 

0 •• 

• 1/ 


9 



where the sequence (r',..., r") of length m — 1 is the sequence (1,2,..., m) with the term 
r excluded. Clearly, det Br = det^,-- Hence 


det Ar = iX 


m m 

Yf ^ ^ (Tp (7 f- 


p=i 
\ P¥=r 


(3.9) 


p=i 
P¥=r / 


and it follows that 

goiX^A) = - 


2^p=i ^p 


L^p=l ^P 


Yr eEi Yp 


Ep=l 


1 ^ r ^ m. 


(3.10) 

(3.11) 


p 


Substituting Eqnations (3.10) and (3.11) into Equation (2.12a), we obtain expressions for 
each in terms of the Eourier transforms of every initial datum Qq^-) and of every solution 


1 r 




/ Er=i qU^) qU^) Ep=i ^p - ET=i € (^ 


+ 


-A \ \r^m 


^r 7 ^ -p ( -A 


ldD+ 


27r 


E m 

p=i ^p 


r gixA/o 

ldD+ 27r 


'Ep=,qH^A 


E m 

p=i ^p 

fi^A) Er=i^p-^.Er=i«^f^;^ 


dA 


E m 

p=l ^p 


E m 

p=i ^p 


dA. 


(3.12) 


Note that we have separated the second integral of (2.12a) into two integrals, which is 
justified provided they converge. However, for A —)• oo from within the closnre of D~^ 
(Z1+), Re(iA) —)> —oo so by Jordan’s Lemma the final integral of Equation (3.12) not only 
converges but evaluates to 0. Hence, we obtain an explicit integral representation of the 
solution. 


dA 

in terms of only the initial data. 

Remark 3.1. In this example we were able to set R = d. This will not generally occur even 
for the semi-infinite rods when the network also contains hnite rods (see, for example, 0 
Proposition 3]). Indeed, i? = 0 is possible here because detAl is monomial which occurs 
for no configurations of rods except the one considered in this section. 




(3.13) 
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0 



Figure 4; m parallel finite rods. 


4 m parallel finite rods 

In this section we consider the case of m G N parallel finite rods of differing thermal 
diffusivity extending between a pair of vertices. Note that we have chosen the rods to all 
be oriented one direction. This is purely for notational convenience, the parameterization 
of the rod from 0 to or from Lj. to 0 is not important. The graph is as shown in 
Figure]^ One of the two vertices is the terminus of every edge, and every edge emanates 
from the other vertex. As \E\ = m, we find it notationally convenient to identify the set 
of edges E with the set {1,2,..., m}. 

The continuity interface conditions allow us to define 


fi(o(A^^) := go (A^,t) , 
/lo(A^t) := ho (A^,t) , 


(4.1a) 

(4.1b) 


for r = l,2,...,m. Similarly, the continuity of flux interface conditions imply 


m 



(4.2a) 


m 



(4.2b) 


r=l 


Using the notation (4.1) the global relation (2.8a) reduces to 


(T,.^Ac/o(A^ t) + c^^5^^(A^ t) 


— 6 


[ariXho{X‘^ ,t) + alh\{X^,t)) 



(4.3a) 


Applying the transformation A i—?■ —A which leaves the spectral functions gj{X‘^,t) and 
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invariant we obtain 

iXL 

-(Jri\go{\^,t) + algi{\^, t) — [—ari\ho{\^,t) + alh\{\^,t)) 


= % 


zh. 

{T y 


Ur 


(4.3b) 


This provides a system of 2m equations valid on C in 2m + 2 unknowns. System (4.2) 


provides two further equations in a subset of the same unknowns. We express this linear 
system as 

AX = Y, 

where 

^ = (5o(A^ t), /^o(A^ t), (A^ t),..., a^g^(A^,t),(^fhl(\^, t},..., , 

Y = (0, yi(A),..., y^(A), yi(-A),..., Wx(a), , 


and 


(4.4) 


>4 = 


/ ° 

0 

1 

1 

... 1 

0 

0 

0 

\ 

0 

0 

0 

0 

... 0 

1 

1 

1 



— iXLi 




— iXL-i 




iXai 

—iXuie "^1 

1 

0 

... 0 

—e '^1 

0 

0 


iXa2 

-i\L2 

—iXa2e "^2 

0 

1 

... 0 

0 

-i\L2 

—e '^2 

0 


iXcTm 

— iXLrn. 






-iXL-m 


—iXo^e 

0 

0 

... 1 

0 

0 

. . . —g Om 



iXL-^ 




iXL-^^ 




—iXai 

iXaie '^1 

1 

0 

... 0 

—e '^1 

0 

0 


-iXa2 

iXa2e '^2 

0 

1 

... 0 

0 

i>'L2 

—e '^2 

0 



iXLm 

iXOmd 

0 

0 

... 1 

0 

0 

iXL-m 

. . . —g Om 

/ 


Solving the system via Cramer’s rule, we obtain 


(4.5) 


5o(A^^) = 


AC -BD 

( B 

A 


/lo(A^t) = 


— det 


Sr{X) 


cr, 


5^(A^^) = 


A 
B 

iXLr —iXLr 

C^(A) yr(A)e^?^~-y,.(-A)e^^ 


BC -AD 
X{A‘^-By 
C \ 

D 


'iXTj’P — 
g 0‘'p _g O'j' 


det 


h[{xyt) = 


a‘^{A^-B‘^) 

( B A D 

A B C 

y^(A)-y,(-A) 

\"^ 5 ,.(A) ^'■Sr(A) 

\ ' ^ O f' _g Of 


(4.6) 


a^{A^-B^) 
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where 


m ^ 


B = Y^a 


Cp(\) 

Pc ’ 




C = 


Yp{\)-Yp{-\) 


i\I/p — i\L-^ 

p=i g — g “"p 


p=i 


Yp{X)e '^p — Yp{—X)e ‘^p 


i\L'n — i\L'n 


g (Tp - g CTp 


for 


Sp{X) = sin(ALp/cTp) 


Cp{X) = cos(ALp/(Tp). 


Equation (4.6) for the t-transformed boundary and interface values depends upon the 
Fourier transform of the solution through Expression (4.4) for Yp. In order for 

Equation (2.12b) to represent an effective integral representation for each solution we 
must remove this dependence. 

The product of a meromorphic function with a holomorphic function must include as 
zeros, every zero of the meromorphic function. As 



2 


(A2 



(4.7) 


is an exponential polynomial (indeed an exponential sum), it is possible to obtain bounds 
on its zeros. Indeed, by m Theorem 3], the zeros of the above product are confined 
within a finite-width horizontal strip. Hence, for sufficiently large R > 0, there are no 
zeros of within U D~^. Indeed, a simple calculation reveals 

(m — 1) log2-|-log(8-|-m(m — 1)) 

R =-^—r- 


is sufficient. A standard asymptotic analysis (see, for example, 0) shows that the terms 
involving qP(-]t) decay as A —)> oo from within the appropriate domains D^, and Jordan’s 
Lemma establishes that these terms make no contribution to the integral representation. 
We have established the effective integral representation for the solution given by 


equation (2.12b) with spectral data given by equations (4.6) where 


y;(A) = go ( ^ ) , 1 ^ r ^ m. 


5 m finite rods connected at a single point 

In this section we consider the case of m G N hnite rods of differing thermal diffusivity 
joined at a single vertex. The graph is as shown in Figure Of the m -|- 1 vertices m are 
the terminus of a single edge and all m edges emanate from the other vertex. As \E\ = m, 
we find it notationally convenient to identify the set of edges E with the set {1,2,..., m}. 
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Figure 5: m finite rods joined at a single vertex 


We enumerate the vertices so that each edge r emanates from vertex 0 and terminates at 
vertex r. 

We use the condition of continuity across the interface to establish 


5^o(A^^) := = / e^''^6'o(s)ds 

Jo 


(5.1) 


for r = 1, 2,..., m, A G C. 

Recall the Robin boundary condition at for each of the terminal vertices with only 


one edge, (1.9). Taking the time transform of (1.9) we have 




(5.2) 


In order to reduce the size of the linear system that will eventually have to be solved, it is 
helpful to use the Robin condition to express either the 0**^ order or the 1®* order boundary 
value at Lr in terms of the other. Suppose that precisely 1 ^ ruAr ^ m of the boundary 
conditions are in fact Neumann conditions (that is the corresponding /3q =0). Moreover 
suppose, without loss of generality, that these Neumann conditions are the boundary 
conditions at vertices 1, 2,..., mAr, that /3[ = 1 for each 1 ^ r ^ rriN, and that (3^ = 1 


for each mAr + l^r^m. We use the boundary conditions (5.2) and Equation (5.1) to 


rewrite the global relation (2.8a) as 


ariXgoiX'^,t) + cr^giiX^,t) — e i\arhQ{X^,t) = 
CTj. J y (Tj. 


9o 


+ e 


alfr{X^]t), 


(5.3a) 
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for 1 ^ r ^ mAT, A G C and 

ariXgoiX^, t) + a‘^gl{X^,t) — e~^ {—iXPl/ur + l)o-lh\{)?,t) = 

or mpf +1 ^ r ^ m and A G C. Another set of m global relations is obtained by evaluating 


equations (5.3) under the map A i—)■ — A, and observing that the spectral functions gj, hj 
and the t-transformed data fr are invariant under this map. 


Together with a t-transform of the interface condition (1.8) the 2m global relations 


form a system of linear equations in the 2m + 1 spectral functions 


9o(A^^), 


t > 0, 

ff^(A^^), 

1 ^ r ^ m, 

t > 0, 

Kix^t), 

1 ^ r ^ mN, 

t > 0, 

hl((A2,t), 

mjv + 1 ^ r ^ m, 

t>0. 


The system can be expressed as 
where 


AX = Y, 


X = {go, (rfgl, ■ • •, crl,g^, crihl,..., al^h^Y, 

Y = {<d,Yi{k),..., Ym{k), Yi{-k),Ym{-k))'^ , 

(i;) “ (^’^) +e ‘’r" alfr{\-,t) ^ m^, 


Yr{\) = 


Qo 




— iXLf 

q' ( —; t) + iXe arfr{X', t) mjv + l^T^m, 


0 li xm 0i xm 

A = I ■To(A) Imxm A\{ A) 

Alo(—A) 

Imxm Ai{\) 
A (A) = (iAiTi, ... ,iAcrm)'^, 

and 

Ali(A) = 


/ iXL^ 

' iXe "^1 


0 


iXe 

0 




— lie 


(5.4a) 

(5.4b) 

(5.4c) 

(5.4d) 

(5.4e) 

(5.4f) 


0 

0 


-a/3r 1A 

-— 11 e y 


(5.4g) 
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The boundary conditions (5.2) give the remaining spectral functions in terms of those 
appearing in the vector X. 

Solving this system and substituting the solution into equation (2.12b) provides an 
integral representation of the solution in terms of the initial and boundary data, and the 
Fourier transform of the solution q^{-] t). It remains to show that it is possible to remove 
the dependence on the Fourier transform of the solution. 

When the system (5.4) is solved using Cramer’s rule, each entry Xj of X is expressed 
as the ratio of two determinants, each of which is an analytic function of A. The denom¬ 
inator is an exponential polynomial in which every exponent is purely imaginary. Hence 
the zeros of the denominator lie asymptotically within a horizontal logarithmic strip m 
Theorem 6]. Moreover, Xj is holomorphic except at the zeros of the denominator, so 


choosing R sufficiently large ensures that Xj is holomorphic on U DJ^. 

As the growing (respectively, bounded) entries of A lie in the same rows as the growing 
(respectively, bounded) entries of Y involving but each such entry in Y grows 

(respectively, decays) like multiplied by the corresponding entry in A, each integrand 
involving g^(-; t) decays as A —)• oo from within the relevant sector D^. Hence, by Jordan’s 
Lemma, the terms involving make no contribution to the solution representation. 


Remark 5.1. Using a similar Jordan’s Lemma argument, it can be justihed that (/J(A^,T) 
and /i^(A^,T) may be used in place of g''j{\^,t) and hj(A^,t) (see, for example, m page 8]). 
Evaluating these f-transforms at final time T only is helpful for numerical calculations. 


5.1 Example 

We study the problem of three hnite rods emanating from a single vertex. The rods have 
lengths 

Li = l, L2 = 1, L3 = 2. (5.5) 

Dirichlet conditions are prescribed at the terminal vertices 

q^{l,t) =sm{t), g2(l,t)=0, q^{2,t)=0, (5.6) 


so that the system is driven by an oscillating temperature profile at the terminal vertex of 
rod 1. The physical properties of the rods are such that their scaled thermal diffusivities 
are 


af = 4, a^ = 9, fj| = 1. 


(5.7) 


For this system, equation (5.4) evaluates to 


/ 0 

1 

1 

1 

0 

0 

0 \ 


/ go \ 

2iX 

1 

0 

0 

_g*A/2 

0 

0 



3iA 

0 

1 

0 

0 

_giA/3 

0 



iX 

0 

0 

1 

0 

0 

_g2iA 


gl 

-2iX 

1 

0 

0 

_e-iX/2 

0 

0 


Ah\ 

-3iX 

0 

1 

0 

0 


0 


9/if 

^ —iX 

0 

0 

1 

0 

0 



\hl) 
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27T < t < Stt 
H Rod 1 
C Rod 2 
□ Rod 3 


t 


7 8 9 



27T < t < 37T 
Q Part of Rod 1 
LI Rod 2 
Q Rod 3 


Figure 6: Three finite rods with 27r ^ t ^ Stt. The heat distributions q^{x,t) (orange), 
q^{x,t) (blue), and q^{x,t) (green) are shown. The code used to produce this plot was 
written by R. J. Buckingham and D. A. Smith. 


The determinant of this system is 


, , , , 7A\ /11A\ /17A\ 

- 4^A ( 2 cos ( — j + cos ( j - 3 cos ( j 


(5.9) 


which has no nonreal zeros so, for any choice of R > 0, the terms involving q^(-,t) make 
no contribution to the solution representation. Therefore, for the purposes of numerical 


work, it is sufficient to solve the system (|5.8|) with 

rT 


Y = 


^0,2iA J sin(s) dA, 0,0, —2zA J sin(s) dA, 0, 0^ (5.10) 


and substitute the resulting formulae for g, h into equation (2.9a). 


Using a Mathematica code written by R. J. Buckingham and D. A. Smith, we obtain 
the plots of q^ (orange), (blue), and q^ (green), for x £ [0,Lr] and t G [0,27r] as shown 
in Figures |5 .1 1 and [b} An animation is shown in Figure [^showing all q^, for 0 ^ f ^ Svr. 
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t= 0 
I Rod 1 
I Rod 2 
I Rod 3 


Figure 7: Please see the ancillary file “threerodsplotfromO-oneshot.swf”. Three finite rods 
with 0 ^ t ^ Stt. The heat distributions q^{x,t) (orange), q^{x,t) (blue), and q^{x,t) 
(green) are shown. The code used to produce this plot was written by R. J. Buckingham 
and D. A. Smith. 



Figure 8: Three finite rods q^{x,t), q‘^{x,t) and q^{x,t) with 0 ^ t ^ Svr. The code used 
to produce this plot was written by R. J. Buckingham and D. A. Smith. 
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6 Concluding remarks 


Remark 6.1. In this work we have provided a general framework implementing the Fokas 
method for solving interface problems for the heat equation on a network. We have demon¬ 
strated the full argument in a number of cases, but it should be noted that the method is 
generally applicable to any such problem. A general argument for the effectiveness of the 
generalized spectral Dirichlet-to-Neumann map, including the Jordan’s Lemma argument 
that terms involving do not contribute to the resulting solution representation, 

was omitted in this paper but will appear elsewhere. 

Remark 6.2. The examples presented in this paper were selected because their general¬ 
ized spectral Dirichlet-to-Neumann maps can be formulated as linear systems which are 
relatively easy to solve by hand. It should be noted that this is not always the case. In 
general, for a graph with mj finite edges and rrii infinite edges, it will always be possible to 
express the generalized spectral Dirichlet-to-Neumann map as a linear system of dimen¬ 
sion 4m/ -|- 2mi but it may not be practical to solve this system by hand. For example, 
the case of m finite rods joined end-to-end in a line (see [5] for the cases m = 2 and m = 3, 
see [2] for arbitrary m with homogeneous Dirichlet conditions) has a generalized spectral 
Dirichlet-to-Neumann map expressible as a pentadiagonal linear system. It is the authors’ 
opinion that a computer code would be the most appropriate tool to obtain the solution 
representation for arbitrarily complex networks. 

Remark 6.3. It has been established that the Fokas method may be used to solve all well- 
posed initial-boundary value problems for a linear constant-coefficient evolution equation 
on a finite interval or half-line miiniiii]. A number of other interface problems have 
already been studied using the Fokas method for the heat equation 0 El uni E], the 
linear Schrodinger equation nznn, and a third-order dispersive equation [6]. Given the 
results above, it is reasonable to expect that the method could be used to solve problems 
for higher order PDF on a general network. 

In the present work, we made use of a particular physical problem to determine ap¬ 
propriate interface conditions but such a physical model may not always be available. 
However, as the Fokas method can be used to determine well-posedness criteria for half¬ 
line and finite interval problems dSlEniEIlES], it should be possible to determine interface 
and boundary conditions that specify well-posed problems for higher order equations on 
networks. 


Acknowledgements 

This work was generously supported by the National Science Foundation under grant 
number NSF-DGE-0718124 (N.E.S.). The authors wish to thank Bernard Deconinck for 
helpful discussions. Any opinions, findings, and conclusions or recommendations expressed 
in this material are those of the authors and do not necessarily reflect the views of the 
funding sources. 


19 


References 

[1] M.J. Ablowitz and A.S. Fokas, Complex variables: Introduction and applications, 
second ed., Cambridge Texts in Applied Mathematics, Cambridge University Press, 
Cambridge, 2003. MR 1989049 (20041:30001) 

[2] M. Asvestas, A.G. Sifalakis, E.P. Papadopoulou, and Y.G. Saridakis, Fokas method 
for a multi-domain linear reaction-diffusion equation with discontinuous diffusivity, 
Journal of Physics: Conference Series 490 (2014), no. 1, 012143. 

[3] A. Bejan, Constructal-theory network of conducting paths for cooling a heat generating 
volume, International Journal of Heat and Mass Transfer 40 (1997), no. 4, 799-816. 

[4] H.S. Carslaw and J.C. Jaeger, Conduction of heat in solids, 2nd ed., Oxford University 
Press, New York, 1959. 

[5] B. Deconinck, B. Pelloni, and N.E. Sheils, Non-steady state heat conduction in com¬ 
posite walls, Proc. R. Soc. A 470 (2014), no. 2165, 22. 

[6] B. Deconinck, N.E. Sheils, and D.A. Smith, Linear KdV with an interface. In prepa¬ 
ration (2015). 

[7] B. Deconinck, T. Trogdon, and V. Vasan, The method of Fokas for solving linear 
partial differential equations, SIAM Rev. 56 (2014), no. 1, 159-186. 

[8] A. S. Eokas, On the integrability of linear and nonlinear partial differential equations, 
J. Math. Phys. 41 (2000), 4188-4237. 

[9] _, A unified approach to boundary value problems, CBMS-SIAM, 2008. 

[10] A. S. Eokas and B. Pelloni, Two-point boundary value problems for linear evolution 
equations. Math. Proc. Cambridge Philos. Soc. 131 (2001), 521-543. 

[11] J. Eourier, Analytical theory of heat, Dover Publications Inc., New York, 1955. MR 
0075136 (17,698f) 

[12] J. Kevorkian, Partial differential equations, second ed.. Texts in Applied Mathemat¬ 
ics, vol. 35, Springer-Verlag, New York, 2000. MR 1728947 (2000i:35001) 

[13] R. E. Danger, The zeros of exponential sums and integrals, Bull. Amer. Math. Soc. 
37 (1931), 213-239. 

[14] D. Mantzavinos, M.G. Papadomanolaki, Y.G. Saridakis, and A.G. Sifalakis, Fokas 
transform method for a brain tumor invasion model with heterogeneous diffusion in 
dimensions. Applied Numerical Mathematics (2014), -. 

[15] B. Pelloni, Well-posed boundary value problems for linear evolution equations on a 
finite interval. Math. Proc. Cambridge Philos. Soc. 136 (2004), 361-382. 

[16] N.E. Sheils and B. Deconinck, Heat conduction on the ring: Interface problems with 
periodic boundary conditions, Appl. Math. Lett. 37 (2014), 107-111. MR 3231736 


20 



[17] _, Interface problems for dispersive equations, Studies in Applied Mathematics 

134 (2015), no. 3, 253-275. 

[18] _, The time-dependent Schrodinger equation with piecewise constant potentials, 

In preparation (2015). 

[19] K. Shivakumar and S. Narasimhan, A robust and efficient NLP formulation using 
graph theoretic principles for synthesis of heat exchanger networks. Computers and 
Chemical Engineering 26 (2002), no. 11, 1517 -1532. 

[20] D. A. Smith, Spectral theory of ordinary and partial linear differential operators on 
finite intervals, Phd, University of Reading, 2011. 

[21] _, Well-posed two-point initial-boundary value problems with arbitrary boundary 

conditions. Math. Proc. Cambridge Philos. Soc. 152 (2012), 473-496. 

[22] _, Well-posedness and conditioning of 3rd and higher order two-point initial¬ 

boundary value problems, 2013. 

[23] T.F. Yee and I.E. Grossmann, Simultaneous optimization models for heat integration- 
II. heat exchanger network synthesis. Computers and Chemical Engineering 14 
(1990), no. 10, 1165- 1184. 


21 



